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It is essential to have these books; the course is based on them and will 
not make sense without them. 
Conventions 
Before working through this correspondence text make sure you have read 
A Guide to the Linear Mathematics Course. Of the typographical conven- 
tions given in the Guide the following are the most important. 
"The set books are referred to as: 

К for Ал Introduction to Linear Analysis 

N for Linear Algebra and Matrix Theory 


AII starred items in the summaries are examinable. 


References to the Open University Mathematics Foundation Course Units 
(The Open University Press, 1971) take the form Unit M100 3, Operations 
and Morphisms. 


Note 


Please note that this text is not based on the set books for the course. 
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31.0 INTRODUCTION 


In Unit 23, The Wave Equation, we used two distinct methods to solve the 
wave equation 


Pu Ou 
Ox! a? Or 


The first method, that of separation of variables, is appropriate when the 
variable x is restricted to a finite interval, say [0, л], and there are boundary 
conditions on u(x, t) when x is at either end of this interval. (Let us note 
that independent of interval and boundary conditions, it is the form of an 
equation which allows us to separate it into ordinary differential equations 
for each variable separately. But it is the interval and the boundary con- 
ditions which determine whether we can find a useful solution thereby. 
Separation of variables gives the solution in the form of a series of 
functions.) 


The second method, D'Alembert's solution, is appropriate when хє А 
and gives the solution in the form 


u(x, t) = f(x + at) + g(x — at), 
where f and g are functions of one variable. 


When we try to solve other partial differential equations, for example the 
equation for heat conduction, which is considered in Unit 32, The Heat 
Conduction Equation, we find that the separation of variables method is 
Still useful. However, D'Alembert's solution is not, being too closely 
dependent on the special characteristics of the wave equation. We have 
already seen a method which can be used in some of these cases, viz., the 
Laplace transform. In this unit we shall discuss another transform method, 
one which js particularly suited to situations involving waves and other 
phenomena concerned with energy transfer. This is the method of Fourier 
transforms. It has other applications besides partial differential equations, 
for example in probability theory, quantum mechanics, and the theory of 
infinite groups, but the only application we consider in this course is to 
partial differential equations. 


The basis of the method is to think of the infinite interval (— co, со) asa 
limiting case of a finite interval [—p, р]. Before trying to solve any partial 
differential equations at all, we consider the simpler problem of applying 
this idea to an ordinary Fourier series. As we shall see, in the limit for 
large p, the Fourier series becomes an integral. 
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311 THE FOURIER INVERSION FORMULAS 
31.1.1 The Fourier Cosine Transform 


Let us carry out the limiting process outlined in the Introduction, for the 
Fourier series of a function f which satisfies certain restrictions imposed to 
simplify the calculation. (Some of them will be removed later, after we 
have established the basic formula.) The restrictions are 


(i) f has domain and codomain А 
(ii) / is continuous and has a derived function f’ which is piecewise 
continuous (page K330). Further, f" has a derived function f^ which 
is piecewise continuous. That is, f is piecewise smooth. 
(ii) fiseven:f(—x) = f(4-x) for all x e R. А 
(iv) There is a positive геа] number ро such that f(x) = 0 for all |х| > ро; 
when a function has this property we say that it is finitary*. 


From Unit 22, Fourier Series (Section 9-6 of K) we know how to represent 
such a function fas a Fourier series in any finite interval [—p, p], where р 
is any positive number. The formula is 
1 X 2nx 
Го) = 5 40 + 4 cos fasse (хєе[-рр}, () 


where 


Р 
«= 1] ficos a 
P l-r Р 


= : [ F(x) cos = dx (2) 


The series (1) contains no sine terms, because / is even, and converges 
pointwise to f(x) for x € [— p, p) because of condition (ii). For x outside 
[— р, р}, however, the series may not converge to f(x); indeed, the sum of 
the series defines a periodic function of x with period 2p, whereas in 
general f is not periodic. 


sum of series 


-àp -p | P 3p x 


Clearly, the larger the value of p, the larger the interval over which we 
get a Fourier approximation. To get a representation that works for all х, 
we take the limit for large р of Formula (1). Then, whatever x is chosen, 
p will eventually be larger than |x| and so x will lie inside [—р, p]; thus 
the formula holds, in the limit, for that value of x. This argument gives 


a ffl X 2i 
10) = lim (в + а e084 a cos Z +) (хє А) 
р 


prem 

(3) 
where. "p co" means the same as “р large." (The justification of 
Equation 3 is considered in Exercise 3 of this sub-section.) 


To calculate the limit (3) we need to know how the а, depend on p for 
large P. This can be found from Equation (2), using the fact that f is 
finitary. For then there is a number Po Such that 


f(x) =0 for |x| > py. 


* Some authors use the word finitary to cover (i), (ii) and (iv). 
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Then if p > ps; we have f(x) =0 for |x] > p, so that in this case, Equa- 
tion (2) reduces to 


2[° & 
а= zf f(x) cos ic dx. 
P Jo Р 
Thus we may write 


2. [Кё 
аг= =F, (©) if p> 4 
= Сүр P> Po (4) 
where 
EI 
ЕДУ) = f Го) cos (ух) dx — (yeR) 6) 
This is an important definition. We call the function F, the Fourier cosine 
transform of f. 
E 
Foy — | S(x) cos (ух) ах — (yeR) 
° 


Using the formula (4) in (3) we can write the limit we want to calculate 
in the form 


2[1 2; 2 
P(x) = lima баю +F; 6) cos + F, (=) cos 27 y E (6) 
pom P \2 р, Р Р р 
When p becomes large, the spacing between the successive points > 
2 
et +++ in the domain of the function F, becomes very small; these points 


cover the positive real line more and more densely. This suggests replacing 
the sum in (6) by an integral. In Unit M100 9, Integration I, we defined the 
definite integral as follows 


LI 
| g= lim A(g(a) + g(a + h) + + gla + (n — 1)4)) 
a noe 
where Л = 2-е It is reasonable to suppose that, for suitably well- 
n 
behaved functions g this formula also holds in the limit for large 5, i.e. that 
а 
| g= lim A(g(a) + g(a + h) + g(a + 2h) +- se) (7) 
а hoo 


This is proved in Appendix | for functions satisfying conditions (i) to (iv). 
We can use (7) to simplify (6); we set 


a=0 
hae 
р 


and 


g(kh) = F, & ex (73) є 


80) = F(y)cos(xy). ` 
Then (6) becomes 


f(x) = lim Zaro + Fe(h)cos(hx) + F.(2h)cos(2hx) + ---) 
h0 
En h(4g(0) + g(h) + g(2h) +) 
T a0 
A EN 
=: m (gO) + g(1)- ---) = lii 380) 


2 (= 
== +0. 
HIE 
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Thus the limiting form of the Fourier cosine series for f is 
2 (= 
Јо) = Ч) Е.(у)соѕ(ху) dy (8) 
о 


A formula to recover a function from its transform is known as an inversion 
formula. In the case of Fourier cosine transforms, Equation (8) is the 
inversion formula. It tells us that the Fourier cosine transform of the 


"n LM " s 
Fourier cosine transform of f is just 5 times f itself. Equation (5) shows us 
how to calculate the Fourier cosine transform of a given /; Equation (8) 
shows us how to recover f if we are given its Fourier cosine transform. 
Example i 
The function f defined by 


р if jx <1 
fa-[. ES) gew 


satisfies conditions (i) to (iv) (Check this.) 
lts Fourier cosine transform is F,, where 


ЕДу) = [^ = x)cos(xy) dx 
o 


0S " " 
ы LY (integration by parts). 


The inversion formula then gives us the complicated integral 


1=|x| if |x| <1 


2 (°l—cosy 
=| — 7 cos(xy) dy = = 
f 2 соо) у sso „л i ee 


л 


Exercises 


1. Follow through an argument analogous to the one given in the text 
above to obtain a formula giving f in terms of its Fourier sine transform 
Е, when f is odd: f(—x) = —f(x). Е, is defined by 


FQ) = [чабу ах (eR). 


2. Find the Fourier cosine transform of the function f defined by 


cosx if Ix «5 
Је) = (хє) 
о if |х|>= 
2 


and write down the integral obtained from the Fourier inversion 
formula. 


fix) 
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3. Using the definition of a limit, prove that Equation (3) follows from. 
Equation (1). 
(Hint: consider an arbitrary fixed value of x, and take the limit for 
large p.) 


Solutions 


1. The changes necessary are: 
Condition (iii) is replaced by 


Gii)’ f is odd 
whence (1) becomes 
21 
Јо) = b sin p sin 4--- (pep p pp 
Р Р 
and (3) becomes 


2 
Уо) = lim Q sin == 4 p, sin 4. ) (кєк) 
p^o Р р 


In the rest of the argument replace 
ао by 0; a, by b, if k #0; cos by sin; and F, by F,, 


throughout. The final result is 


го) =$ [онтоо as 


x2 
2. Еу) = | cos x cos (xy) dx 
о 


= n [1 cos(1 + y)x + 3cos(1 — y)x] dx 
° (Section III.5.1 of TI) 
in x sin(l—y»x| , 
geet mer е on 
Sin dz(l +y) , sindz(1 — y) 
20 + y) 20 — y) 
cosiny  cos(—imy) 
T+A 20-5 
_ cos јлу 


- T. 


1-y 


For y = +1, a separate calculation gives 
к 
1) =-. ' 
ЕДЕ) = 


The inversion formula gives 


x 


cosx if |x| *3 


л 


2 L cos(1xy)cos(xy) РА 
— үе эса 
$ á о df ЕЕЕ 
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3. Equation (1) tells us that 
Јо) = Сір, х), it xe[-p, р] (9) 
where 
C(p, x) = ; do + cos (=) " 
We wish to show that 
fG)-limCpx) (eR). 
27 


i.e., for any x e R and for any positive number c there exists 


а real number № such that : 


|70) – С(р, x)| < £ for all p > N. 


We can achieve this by choosing № to be |x|, for Equation (9) 
tells us that 


ГО) — Cip, x) = 0 for all р> |х|, 


and 0 is certainly less than c. Since this argument works for 
any real x, we have proved (3). 


31.1.2 The Connection with Laplace Transforms 


In the preceding sub-section we defined, for suitable types of function f, 
the Fourier cosine and sine transforms F, and F,: 


Fo- [оно oem a) 


and 
F0)= ferens oem o 


Although we stated that f must be even in the first case and odd in the 
second, the actual integrals only involve f(x) for x 2 0, and so we can use 
these two expressions to define Fourier sine and cosine transforms even 
if f only has domain [0, со) instead of (— co, со)—ргоуійей, of course, 
that the integrals converge. Of course, if f has domain [0, оо), it cannot 
be said to be either even or odd. 


There is a similarity between the above two formulas and the formula 
defining the Laplace transform of f (Unit 29, Laplace Transforms): 


ШЛО) Е = 


Tn ай! three cases f(x) is multiplied Ьу a particular expression involving y 
and x and the product is integrated. Since in each particular case the 
result depends on y only, it defines a new fünction (of y); such a new 
function is called an integral transform of f. Integral transforms provide a 
useful technique for solving various types of mathematical problem, as 
we have already seen in Unit 29. There are many other types of integral 
transform, for example the Mellin transform defined by 


[о x7 de; 
A 
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or the Stieltjes transform 


e f(x) di, 
o x+y 


but we shall be concerned only with Fourier and Laplace transforms in 
this course, In all cases, we are dealing with a linear transformation on 
some relevant vector space of functions. The utility of these integral 
transforms depends upon the form of the differential equation to be 
solved. 


Exercise 


If К is а continuous function from R x R to R show that 
b 
A fr— | Kix, yf(y)dy (єў) 


is a linear transformation on V, the relevant real vector space of functions. 
(Note that K might be a Green's Function, Laplace transform, Stieltjes 
transform, or almost anything.) 
Solution 

We show that 


Alaf + fg) = «А(/) + BAC). 
Alaf + Bg) = [к (х, af + fig) dy 


b 
= | Ко, yas) + Be dy 
b 
E [e аро) dy + || Ко, уйе) dy 


=a fre »/0)dy +В [к DEO) dy 
= «А(/) + PA(g). 


As you may have suspected from the similarity of the formulas defining 
them, a connection between the Fourier cosine and sine transforms on the 
one hand and the Laplace transform on the other may be found. This is 
because the sine, cosine and exponential functions are themselves con- 
nected by Euler's formula, which you met in Unit M100 29, Complex 
Numbers IT: 


cos Ô + i sin @ = e? (4) 


This formula suggests that we consider the complex linear combination 


БУ) iE) = IN [cos(xy) + їзїп (ху)] dx 


= Јо dx i 
0 
-f£[/K-iy) 6) 


This formula provides the connection between Fourier and Laplace trans- 


forms. There are several points to be careful about here. 


(i) The formula (5) only works if f has domain [0, co) and all the 
integrals converge. 

(i) Тһе integrand is now a special sort of complex function; its codo- 
main is C, the set of complex numbers, and its domain is R. We 


п 
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define the integral of such a complex function by treating the real 
and imaginary parts separately: if g = h + ik, where h and k both 
have R as domain and codomain, then 


| р = | "à +) = f +i [s (6) 


For example, the real and imaginary parts of the function 


8:0 ———— ехр(іб) (8 є [a, 5]) 


£:0:.— —5 cos 0 + isin 6 
are respectively 


6.—. 5 cos8 and 0 — —» sin 6 — (0e [a, b]). 
Thus 


b 
| ехр(і0) 40 
means 


b b 
[oi [sino ao 


= (sin b — sin a) + i(— cos b + cos a) 
_ exp(ib) — exp(ia) 
7 Я 


(ii) We originally defined £[/] to have as domain some subset of R 
(sub-section 29.1.3 of Unit 29). Now we are extending the domain 
to include some complex numbers, such as iu with u real, as well. 


From the point of view of this course, a major difference between Fourier 
and Laplace transforms is that to invert Fourier transforms we have a 
general method based on integrals such as 


2 оо) dy 0) 
т Јо 


Whereas to invert a Laplace transform we have had to use the special 
techniques and tables discussed in Unit 29. This difference is not as funda- 
mental as it may look, because an inversion formula for the Laplace trans- 
form does exist; however, we shall not study this formula because it 
involves a type of integration which is outside the scope of this course. 


Example 
Consider the function (note the domain) 
J:x e — (xe[0, oo) 


where a is some positive number. Note that f is not finitary, but even so 
the integrals defining F, and F, converge. 


Working out the integrals that define F., Е, and £[f], we obtain 
Fo) Í e^ cos(xy) dx 
о 


me 
pru 


OER), 


FQ) = | ке sin(xy) dx 
o 


з OER, 
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апа 
ШО enm 
o 
тел. 
“зүп (sao) 
Thus, 
И 1 
ОЛ) = ETET 
iy+a 


“Cit UCET) 
ЗС LR 
= ууа 

50 that Formula (5), 

ЛС) = F.0) + iF O), 

is verified. 

Example 

For the function 

fix——— e (x € [0, со)) 


with a again positive, we have 
f 
CL) =- (e(a, co) 
у-а 


but F, and Р, are not defined, since the relevant integrals diverge. Indeed, 
the condition: 


E 
Í е°* dx converges 
о 


is not satisfied. We can still calculate the real and imaginary parts of 
EK- iy) 


but we would not be justified in applying (5) and calling them F.(y) and 
Fy) respectively. 


Exercises 
ъ If 
J: xı 2ix-e* (хє К), 
calculate 


3 
Í Дх) dx. 
° 


2. (i) Find the Laplace transform of the function 


fix— [73 p i] (x € (0, co). 


(1) Evaluate the real and imaginary parts of 
E-i). | 
Ате we justified in calling these F.(y) and F,(y) respectively? 


(iii) Check by comparing with the formula for the Fourier cosine 
transform of a function agreeing with f on [0, co), which we 
calculated in sub-section 31.1.1. 
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Solutions 


14 


E 


2. () 


Gi) 


(iii) 


ЯЯ 
| (2ix — e?) dx 
о 
= ("cos 3x) + ох — sin 339] dx 
a 


е чт . 
- [i-em 3x]dx + if (2x — sin 3x) dx 
о 


o 


sin 3x]*/? | 2, Cos >" 
= + |e 
| 3 ji DT 3s 


1 (v 3) 
Calculation of C[f (y) ' 


[ro de = [о = e dy 
А 


o 


e"? pe 22) 1 
КЕ -=y Y jl 


y 
et od ot 
y ty y 
Calculation of the real and imaginary parts of €f ( — iy) 
А е?” 1 1 
VIT уу a Cg 
_cosy+isiny i t 
ту y ¥ 
_ 1 22-50) 
У у y 
Application of the formula 7] = Е, + iF, 
We have 


FAy) = [л (x) соз (xy) dx 


= [о — x) cos (xy) dx 
o 


which converges since the end-points are both finite. 
The integral for F,(y) also converges for the same 
reason and so we may conclude that 


1 — cos 
Fy) = =F 

1 siny 
КО = y y 


Comparison with direct calculation of F. 
In sub-section 31.1.1 we calculated the Fourier trans- 
form of the function 
1—|x| if |x| «1 

0 # |> 26 
Although this is not the same function as our f its 
images agree with f(x) for x > 0 and hence the trans- 
form F.(y) is the same for both functions, We obtained 
there, by direct integration, 


—| 


1=cos у 
Еу) = — = 


which agrees with the result obtained above by the 
Laplace transform. 
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31.1.3 Complex Fourier Transforms 


After our digression in the preceding sub-section about functions with 
domain [0, оо) we return to our main preoccupation, functions with 
domain R. We wish to remove the restriction made in sub-section 31.1.1 
that the function f must be either even or odd, while retaining the useful 
symmetry in the formulas relating a function and its transform. More 
precisely, we have different definitions of the Fourier transform for func- 
tions with domain [0, со) and R; they are related by defining, for f with = 
domain [0, co) an even and an odd extension (see pages K350-1); for 
example, the even extension Ф is defined by 

Ф: Ri К 

Фо) =f(x) if x20 

фо) = ф(х) if x«0. 
Suppose, therefore, that f is a function that satisfies (1), (ii) and (iv), but 
is not necessarily odd or суеп. Since we already know how to deal with 
even and odd functions, we start by expressing f with domain R as the 
sum of an even and an odd function, using the method given in Unit 22, 
(Section 9-3 of K). 

f 7 fe fo 
where 

Ж) = 0) + /(—х)) 

Јобх) = 4) — Л). 
We can then calculate two transforms: the cosine transform of fe(x), and 
the sine transform of fo(x). 


FAY) = [o €os (xy) dx 


FQ) = [so sin (ху) dx 


It would be possible to regard this pair of functions (F., F,) as the "*trans- 
form" of f, and we could use the inversion formulas from sub-section 
31.1.1 to recover fe, fo, and hence f when this pair of functions was 
known; but this procedure is cumbersome and lacks the pleasing symmetry 
of the formulas we gave in sub-section 31.1.1 for the cases where fis even 
or odd. What we would like is a useful way of combining F, and Е, into a 
single function in such a way that this new function is related in a sym- 
metrical way to f. 


A clue to the way to do it is provided by the preceding sub-section, where 
we saw that it is profitable to consider the combination F, + iF,. Let us 
try the same thing here. The combination we are interested in is 


Еду)  iF() 
m NT di Í "folx)sin(xy) dx 
5 


o 
-p tfc Јо) -f(-x) 
2 


2 cos(xy) + i 


sina») dx 
o 


= [ (В/О) сокбху) + i sinfx)] ' 
+ MC X) cosy) — їзї(хуўх 
- ihe Ge + 1f(- x)e7 d 
(by Euler's formula) 
- I (even part of [/(х)е'*^*]) dx 


o 


М 
= Al S(xe* dx — (Equation (9-5), page K335). 
5s 
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Once again we have an integral transform, differing somewhat from the 
others we have considered in having range of integration (— co, оо) in 
place of [0, со) and in being complex for real y. We call 


1f 
> | Jf xe" dx 
2j-» 

the Fourier transform of f. 


We can find its inversion formula from those for the sine and cosine 
transforms as follows 


Је) = f) + foo) 


2° 29 re 
=- j| F(y)cos(xy) dy += Í F,(y)sin(xy) dy 
z Јо n Jo r 
2 isy gn ler ws 
=; отар 


(by Euler's formula) 
"ye? (E) = iF(9)) + 3e" 'P[F (y) + iF(y)] dy. 
o 


If we write F(y) = F.(y) + iF,(y), this becomes 


Је) = E NL + de" F(y)) dy 


o 


F(Y) dy 


1 
т 


I9... 
= чау, 
z [ xS F(y) dy. 


e F(y) dy 
о -0 


L] 


Thus for Fourier transforms, as for sine and cosine transforms, the 
inversion formula is very similar to the formula defining the transform 
itself. 


It is usual to make the formulas even more symmetrical, by re-scaling the 
definition of the Fourier transform as 


ad P" эрлан» 
ЈО) = gir dx 


The inversion formula is then 
1 ° T 
= -ixy 
fe) =F Г Je» dy. 


These formulas are the central ones for us. The re-scaling has a role other 
than symmetry. If f belongs to a certain Euclidean space, the Fourier 
transform is length preserving with this scale. We discuss this further in 
Section 31.3. 


Аз a further extension of these ideas, we would like to do away with the 
condition that f be finitary. It turns out that this can be done. A sufficient 
condition for the validity of the above inversion formula is that f be 
continuous and piecewise smooth and that f be absolutely integrable, by 
which we mean that the integral 


[` лө ж 
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is convergent. If f Баз a jump discontinuity at x = a, say, it can be shown 
that the inversion formula converges to 


Га?) fa) 
P 


provided f satisfies certain conditions. (Compare this with Theorem 9-1, 
page K340.) 


Example 

Find the Fourier transform f of f, where 
Јо) =e — (xeR) 

and a is some positive number. 

The Fourier transform is f, where 


ЈО) = 75 Г етеін dx 
л J= 


= xl е7°\*\ cos (xy) dx + xl ет! sin (xy) dx 


Jn 


since 

xı el?! cos (xy) 
is even and 

xı el sin (xy) 


is odd. Using the result of the example in sub-section 31.1.2, we can 
evaluate the integral, obtaining 


2a 
=>. 
Vala + у?) 
Applying the inversion formula we obtain the integral 


ET 


Tix dy 
a 


Exercises 
1, Calculate the Fourier transform f of f, where 
f(x) = e^*! sin bx (x e R). i 
2. Show that if: 
J: R— R 
is even then the image set of fis a subset of R. 
3. |f 
fiR—R 


is odd, describe the image set of f. 
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Solutions 


1. Јо) = |. e? sin bx е” dx 


eas I: е sin bx (cos (xy) + i sin (ху)) dx 


DL [ e7* Pli sin (bx) sin (xy) dx, 


van Jo 


(since the real part of the integrand is odd and 
the imaginary part is even) 


2 M cos (b — y)x 


prit 


— cos (b + y)x] dx 


EN. a Е а ) 
Ji EED 


by the formula in the example of this sub-section. 


2 fo) = Tz [ze (cos (xy) + i sin (xy) dx 


-——— Г I(x) cos (xy) dx 


(by Equation 9-6, page K335, since f(x) sin (xy) is odd). 
Since /(у) is real, the image set of fis a subset of А. 
3. If f is odd, f(y) is imaginary (i.e. its real part is the zero 
function for all y є R). This is because 
X > f(x) cos (xy) 


is odd and hence 
Г f(x) cos (xy) dx = 0. 


The image set of / is therefore a subset of the imaginary 
numbers (it : t € R). 


31.1.4 Properties of Fourier Transforms 


In formulating the definition and inversion formula for Fourier transforms, 
we have been led to define a new type of integral, in which the integrand 
is a complex-valued function (i.e. one whose codomain is C, the set of all 
complex numbers). For such functions we can integrate the real and 
imaginary parts separately, but the manipulation of such integrals is much 
simplified if we do not timidly convert them at once into integrals of real 
functions, but instead use complex numbers throughout. For example, the 
result we obtained in sub-section 31.1.2 * 


s 
f e" qg = 
a 


could have been obtained by assuming that the formula 


р е ах = G [e — e] 


applies even when z is complex. For example, if F is defined by 
F(8) = ё = cos 8 + ising 


а) 
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we have 
F@= 


as we might expect from the formula 
d 
L (өгү шов 
as (e%) = ae 

for real a. We also have 


b 
|; ie? d9 =e — ее 


" 
and the result (1) follows on dividing both sides by i. 

Since the fundamental theorem applies to such complex integrals, the rules 
of integration derived from it apply too. In particular, we may use inte- 
gration by substitution and integration by parts for complex integrals 
just as for real ones. There is one caution, however: substitutions must 


always be done with real-valued functions (see Example 1 below). This 
will ensure that the new integrands have domain R. 


Example 1 (integration by substitution) 
Given that the Fourier transform of f is J, what is the Fourier transform 
of the function 


xr—+ f(x — a) (xe А) 
where a is a real constant? 


We are given 


j ГО) 0 dx 


and we want 
1 fe 
xl Jf (x — a) е? dx, 
„2л -0 


ICA: 
We make the substitution x — z -- a obtaining, since an 1, 


а as Каду dz 2; 
pE 0) 


The integral (2) can be written 


ei [а 
me | fe"? dz 
2n J-» 
= e fly). 
The Fourier transform of 
x f(x — a) 
is therefore 
y —— e fly), 


This is the analogue for Fourier transforms of the second shifting theorem 
for Laplace transforms (page K195), which states that if a is positive and 
if f is defined by 
1-а) if t»a 
f0- fK | if t<a 
then 
£[/] =e "t(g]. 
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Notice, however, that if a had not been real the method would have 
broken down; in that case 2 = x — а would have been complex, so that 
the transformed integrand in (2) would have had a complex domain, and 
wc have not discussed the integration of functions with complex domain. 


Example 2 (integration by parts) 

Given that the Fourier transform of fis f, find that of its derived function 
У". (Assume that the derived function f’ has a Fourier transform and that 
f(x) approaches 0 for large positive and large negative x.) 

The Fourier transform of f’ is f’ where 


= Z Í ж Го) e dx. 
Integration by parts gives 

"(у= ET "rr -4f ye? dx к) 

ғо) Б Кас, Sye (3) 


-® 
where [ф(х)]®„ means, for any function ф with domain R and codomain C, 
lim ф(х) — lim ф(—х), 
хло x^ 
provided both limits exist. Since we are assuming that f(x) approaches 0 
for large positive and large negative x, (3) reduces to 


i fe 
o=- = | Sye? dx 
Ут )-= 
and so the Fourier transform of f” is 
Piye— yo) 
(What commutative diagram may be drawn here?) 


There is another technique of integration which is useful for Fourier 
transforms, called differentiation under the integral sign. We introduce it 
now merely because we have not had occasion to use it before, and not 
because of any logical connection with Fourier integrals or complex- 
valued functions. The rule states that, under certain conditions, the 
differentiation linear transformation with respect to one variable commutes 
with the integration linear transformation with respect to another: 

d fè ins b дф . 

sl (x, y) dx = [ ay 04 (4) 


where ф is a suitable function from R x R to R. 


Example 3 
If ф(х, y) =e”, then the left-hand side of Equation (4) is 


The right-hand side of (4) is 


ъа à 
= = » 
T ET (e?) dx [^ dx 


5 
(by Section IIL.5.2 of TI) 
а 


be” — ae? 
y 


as before. 
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The formula (4) for differentiation under the integral sign may be applied 
> д 
whenever the functions ¢ and 2 are continuous. (See Theorem 1-35 on 
y 
page K670.) If the interval of integration is infinite, then the formula is 
the same: 


d (= * 2 
5 [7 вела -[ оза 


but the conditions for it to be valid are more complicated: it is sufficient 


д 
that ф and z be continuous; that there exist functions ф, and ф, such that 


1óGs, у)| < (16) and < ol) 


% 
ду 65 у) 
for all у; and that 

E dite) de ma [^ dy) ax 


both converge. A proof is given in Theorem 1-41 on page K676, but you 

аге not expected to study it. 

Example 4 

Given that the Fourier transform. of fis f, what is the Fourier transform of 
8: x xf (x) (хє А)? 

(It is assumed, of course, that both f and g have Fourier transforms.) 


The required Fourier transform is 


2) = ж © эде» 


1 


2. [9 ды 
= = ә ede 


2 ы еу 

= gaje dx 
d 

= - 12570) 

-iyo). 


These four rules: the fundamental theorem of calculus, integration by 
substitution, integration by parts, and differentiation under the integral 
sign, provide the basis for the manipulations of complex Fourier trans- 
forms that we shall be performing in the rest of this unit. We have stated 
conditions for their validity, but they are often difficult to verify in practice. 
We may note that there are other conditions, not equivalent to the ones 
we have given. The point is that the conditions are related to the space 
the function belongs to. What is sometimes confusing is that the same 
function may belong to different spaces! We examine Euclidean Spaces 
in Section 31.3. 


Li 


Example 5 


The Fourier transforms of certain functions are particularly useful in 
applications. One that is especially so is the function 


fixe—— exp(-ix?) (хед) 
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associated with Gauss. To find its Fourier transform we evaluate the integral 
Ў) = L Г exp (—$27) e'? ах. (5) 
J22 )Ј-ә 


We can do this by means of the following device. We differentiate both 
sides of (5) with respect to y, using the rule of differentiation under the 
integral sign, which is valid in this case. This gives 


(0) = 


1 (= д 


ue 


exp (—1x?) 5 (е ®)ах ` 


ә Р 
exp ( 41x?) ixe? dx 


==” 2 Ciapc pepe dx 
Vn J - dx 


Integration by parts gives 


üyo- E [-iexp (—4x?)e"™]?_, 


ж 


+ p exp (—42x?) E (e) ax} 
ВА Ы ГА” 
=0-—= exp (—4x7) e"? dx 
Al. 


= -yf(y). 


Thus we have obtained a differential equation for /. Solving it by the 
method described in Unit 4, Differential Equations 1, (page K96) we obtain 


ЈО) = A exp (—4y) (6) 
where A is some real number. 


We can evaluate 4 by using the Fourier inversion formula, which gives 


Ie) = Fe NL 


=A = E exp (—1y?)e7** dy 
=Af(-x) Бу (5) 
= 42 ехр(- х?) by (©) 


Thus we can conclude that A = 


+1. Trying the special value у = 0 in (5) 
gives 


Foe ane азд 
ба |= 


and so, by (6), 4 must be positive. Tt is therefore 


+1, and we conclude 
that the function 


xı exp (—4x?) 


is its own Fourier transform! 
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Exercises 


1. Use the result 


1 fe 
Јо) = [| exp (—$х)е” dx = exp (—3у°) 
/2х J-a 
in Example 5 to evaluate 
ә 
f exp (—4x?) dx. 
2. Ifg and fare related by 


8(х) = /(0х)  (xeR) 


with a some given positive number, how аге & and related? (Hint: 
look at Example 1 again.) 
3. Usethe results of Examples 1 and 5, and Exercise 2, to find the Fourier 
transform of 
J: x exp (—1с(х — 6)?) (хє А) 


where c and b are given real numbers, with c > 0. 
4. Use differentiation under the integral sign to derive from 


х leo 
fiera- E 
0 5 


formulas for 


х х 
(i) | te^" dt and (ii) | Pe dt 
о о 


(Hint: what is H (еу?) 


Solutions 


1. If f(x) = exp (—1x?), we have 


E exp (—4x?) dx = ‚/2л ЈО) = 2n. 
-% 
2. By definition 
w= Г Ладе?" dx, 
NAAR 
Using integration by substitution with x — 2. 


а 


= z 72) by definition. 
a а, 


3. We have 
F(x) = g(x — 5) i 
where 
g : u exp (23e) 
and 
FOREN 
where 


h: v exp (402) 
We must have 
ЈО) = е?% (у) (Example 1) 


23 


LM 31.1.4/31.1.5 


and 
50) = ET (5-) (Example 5 and Exercise 2) 
ve Me 
and so 
4. (i) 
(i) 


31.1.5 Summary of Section 31.1 


In this section we defined the terms 


finitary function (page C6) Se 
Fourier cosine transform (page C7) ы. 
inversion formula (page C8) t. 
Fourier sine transform (page C8) ex 
integral transform (page C10) * 

Fourier transform (page C16) Aj dn 
absolutely integrable function (page C16) *o* 

Notation 


Е. (раве C7) 
F, (раве С8) . 
P4 (page C16) 


"Techniques 
l. Integrate functions from R to C by applying the rule m 
o à b à 
| є-| Са) naif k 
а а a a 
where g = h + ik. 
2. Apply the following techniques of integration to functions from R + 


to C: 
(i) integration by substitution, 
(i) integration by parts, 
(iii) differentiation under the integral sign. 
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31.2 APPLICATION TO DIFFERENTIAL 
EQUATIONS 


31.2.1 The Method 


The application we shall make of Fourier transforms in this course is to 
the solution of partial differential equations such as the one-dimensional 
wave equation (Unit 23) 


Pu 1 8u 

TPI ® 
and the one-dimensional heat conduction equation (Unit 32) 

ðu 1 du 

Ox a? Ot 


when the space variable x is not restricted to a finite interval. In this unit 
we shall apply the method to the wave equation, obtaining D'Alembert's 
solution which was given in sub-section 23.4.1 of Unit 23. 


We shall also show how the method can yield new results, by applying it 
to the heat conduction equation. 


To solve the wave equation by Fourier transforms we start by taking 
Fourier transforms of both sides of Equation (1) with respect to x; that 
is, we write 


ә 9, 
xl aio de! dx 
л /-® 
__% P 1 au 
я. Jo а at 


To clarify the manipulations that follow we introduce some auxiliary 
notation. For each real number t we define the function of one real variable 


Q: x au(x) — (xeR) @) 


(х, De dx (2) 


For example, if t = 2 we have the function 
$3: x ——9 u(x, 2) (хє А). 


The reason for this definition is that we are applying the Fourier transform 
to x only: x and ! are treated differently. Using this notation the left-hand 
side of (2) becomes 


1 f 

= | d. (де? dx 

Уп =% 
the Fourier transform of $,". Now we saw in Example 2 of sub-section 
31.1.4 how to show, using integration by parts, that the Fourier transform 
of à, is 

y —— -0ф0) 
where 
(xe dx 
is the Fourier transform of $,. By applying the same method a second 
time, we see that the transform of ġ,” is 

y —— (ccip'é) = -y 6 
and so the left-hand side of (2) is —y74,(y). 


To simplify the right-hand side of (2) we use a similar idea. We define 
functions, one for each real number y, by 


ууа. фо) (te R) @ 
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So that 
00 = &0)- X | © als, De de. 


The derivative of this function is V, (f), where 


djl ә iam 
ууа) = T Uz E, u(x, t)e ax) 


-— ba es De? dx 


2n J-o 
assuming that the conditions for differentiating under the integral sign 
are satisfied. In a similar way we obtain 


. 1 fe aw 
wO zl P s ne dx, 


Therefore 
4 V, (t) = right-hand side of (2). 


Using our new forms for both the left- and the right-hand side, Equation 
(2) now becomes 

-ay' $9 =H," 
Using (4) we can write this as 

ау) = uy 0 G 


Thus we have reduced the partial differential equation in v, Equation (1), 
to a whole set of ordinary differential equations іп the functions y,, one 
for each real number y. These functions are related to и by 


1 fe 
WO = 810) = Toe [^ a nas 
NAE 

Taking the transform has replaced the double differentiation with respect 
to x by a multiplication by — y?. 
The ordinary differential equations (5) are easily solved. We know from 
Unit 9, Differential Equations II, that the functions 

t sin ayt 
and 

t > cos ayt 


form a basis for the solution space (this is true even though the solution 
space is now a complex vector space). 

Since we are interested in complex-valued functions, it is more convenient 
to use the alternative basis 


{t el t gen) 
and so the general solution of (5) may be written 
P= Ge" Fe" " — (eR) 


where G, and F, are arbitrary constants. For each real number y we may 
choose G, and F, arbitrarily and still obtain a solution of (3). To describe 
the way G, and F, depend on у, let us introduce two new functions 


fiy—F, VER 
&:y-—>G, (ye R) 


The notation prejudges the issue: we shall be interested in the inverse 
Fourier transforms of the two complex-valued functions f and g which 
are at present more or less arbitrary. 
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Me now want to use our knowledge of у, to-find и. Equation (4) tells us 
that 


3,0) = $0), 


so that 


$y) = Ne Je" — (ye R) © 


holds for every t. 


The inverse Fourier transform of 6, is u(x, t), so that 


u(x, 1) = Ж | © 50) dy 


LS E. 
SS Í eTe Ey) + eR] dy by (6) 
м 


We can rewrite this as 


"mc 
u(x, г) al 80) dy 


Jd е nid 
m la e etm fy) dy 0) 
=9(x—at) + f(x + at) (8) 


where g and fare the inverse Fourier transforms of g and f. Since z and f 
are more or less arbitrary, so are g and f—except, of course, that we would 
like them both to be real so as to make u(x, f) real. Equation (8) is just 
D'Alembert's solution of the wave equation, which we studied in sub- 
section 23.4.1 of Unit 23. But as a concomitant of all our work, we also 
have Equation (7), which is often useful. 


Exercises 


Describe briefly the main steps used in the Fourier transform method 
of solution of a partial differential equation. 


2. Using the Fourier transform method obtain a solution for the heat 
conduction equation 
ди 1 ёи 
x? a? at 
where 
u: Rx [0, о) —— R 
ди бз 
and а> 0. (Assume that и and os approach zero for large positive 
ux 
and negative x, and any other existence condition you need.) Give 
your answer in the form of an integral involving an unspecified 
function g, analogous to Equation (7) above. Check your final answer 
by substituting in the above differential equation. 
3. Show how the unspecified function g in the previous solution can be 
determined from a knowledge of the function 
x——+ u(x, 0). 
Solutions 


l. The main steps are: 
() take Fourier transform of the partial differential 
equation; 
(й) by defining appropriate functions of one variable, write 
down an equivalent ordinary differential equation; 
(ii) solve the ordinary differential equation; 
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(iv) generalize to functions the arbitrary constant(s) in this 
solution and use it to obtain the solution of the partial 
differential equation by applying the inverse Fourier 
transform. 


P.D.E. ——— O.D.E. 


| | 


Solution +—— Solution 
The transformed equation is 
1 (= „дши 1 (9 iayl ди 
е i e pA 
Nn J-o Ox’ Рп J-» a^ Qt 


The equation analogous to (5) is (using the same notation) 


(x, i) dx = (x, 0) dx. 


which holds for every real number y. Solving this as an ordinary 
differential equation gives 


00) = фу) = c exp (а?у?) 


and since the arbitrary constant c, may be different for different 
y the general solution of the differential equation for 6, is 


$,0) = &(y) exp (ay?) 
where g is an arbitrary function from A to C. Taking the 
inverse Fourier transform, we obtain the general solution 


"T A [o exp C7 a^^ — ixy) dy, 


For the time being, this cannot be simplified further. 
Check: differentiation under the integral sign gives 


2, 
EE 
1 (=. А 
T f 20) exp (-a* y? t — ixy)(—y?) dy 
Е a 
ôu 
a д 


= 75). 20) exp (—a*y*t — іху)(– а?у?) dy 


so that 


ax? 
as required. 


The previous result, with t = 0, gives 


u(x, 0) 


a 
Ee dy 
а 


Thus u(x, 0) is the inverse Fourier transform of g and hence 
8 is the Fourier transform of x ——+ u(x, 0). The required 
formula is therefore 


z 1 
#0) = WI po 


u(x, буе” dx. 
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31.2.2 Repeated Integrals 


Both the Fourier transformation and its inverse mapping аге linear trans- 


formations (see the exercise on page C11): the two diagrams below are 
commutative. 


(ЕЮ — fag 


GB —— fre 


scalar multiplication 


f ————— у 


| 


7 ———— af 


scalar multiplication 


On the other hand, taking Fourier transforms is not compatible with 
multiplication of functions. 
But is there a binary operation [7] on the space of functions so that 


(fg) 


fxg 


| 


$8 > Jn£s-Uxg 


does commute? If so, finding Г] would enable us to express (f x g)~ in 
terms of f and g above. 
Useful as this would be, finding a binary operation + such that 


C58) f*g 


| 


(ЛӘ —— Jx 


commutes would be even more useful. For we could then exploit our 
knowledge of f and g. For example, from Exercise 2 of sub-section 31.2.1 
we have the solution of the heat conduction equation in the form 


1 ” А 
u(x, t) = —= | #00) exp (а?у? – ixy) dy. 
-e 


Ta 


J: y exp(-ay!) (ує А), 


ux, ) = X]. 80) x (he dy; 


and if the inverse transforms of g and /,, g and f,, are known, and if + 
exists, then 


u(x, 1) = g(x) * (х). 


In this sub-section we shall find the answer to this problem, namely to 
find л in terms of f and g if we know that fi = f x g. That is, we have 


1 = 2 m 
A(x) = =z | TOR) dy а) 
J2n J-0 
as the input datum. But we may replace fin this: 


A(x) = a fe ајот" ау 0) 


A. s 
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The right-hand side is the integral of a function which is itself defined, as 
an integral. Such an expression is usually called a repeated integral. 


An important property of repeated integrals is that we can usually reverse 
the order in which the integrations are done. (This reversal of order has 
analogues in situations we have seen already. For example for an n x m 
matrix 


Y Sa, 


1 с 


since we can calculate the sum of all the elements in the matrix either by 
calculating the л row sums and then adding them, or by calculating the т 
column sums and then adding them.) For finite ranges of integration such 
interchanges 


f | [ Ех, у) 4| dx = [i Бх, у) в) dy о 


are always valid provided only that the integrand is continuous. For 
infinite integrals 


| -f | КСЕ Y Il Fe, 3) Ф) dy в) 


they are valid if F is continuous and if 


[ai (fue. »l 4j d 


converges. (A related result is Theorem 1-40 on page K675.) 


Assuming that these conditions are satisfied by the integrand in the 
repeated integral (2), i.e. by 


SORY)» 
for each x € R, we may use (3)' and obtain 


fo 1 fo " 
ta) = = |” fol | Hoyer dy) de 
Jn J-0" \J2rJ-o 
° 
=— S(2)g(x — z) dz (4) 
NL SE 

The change of variable u — x — z yields an alternative way of writing the 
same result 


= z ig Дх ви) du (5) 


2n J 


It is conventional to use the same variable of integration in (4) and (5). 
Hence (5) becomes 


м) = — [^ fæ- dge) ae © 


уж) 


The binary operation + is defined by (4) or (6). 


x pe e 
Гев) = Tu [ Og - 2 dz = Ж [| _ I0 7 280 de. 


This is an important result called the 
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Convolution Theorem 
If 


h-fxg 
then 


h-f*g. 
Ais called the convolution of f and g. 
Exercise 


In Exercise 2 of sub-section 31.2.1 we arrived at a solution of the heat 
conduction equation: 


ux, ) = = ie EO) exp (~a?y?i)e dy. 


From Exercise 3 of sub-section 31.1.4, the Fourier transform of 


y ——5 exp (-aty?t) 


for each ге К. 
Verify by differentiation that 
u=ged, 


| 3 Ou 1 aw 
(ie g * Q, satisfies э" a) 


Solution 
Using differentiation under the integral sign we obtain 


ди _ oe 1 (z — x)? 
T eau [ s0 e- da ) 


- | &(z) am ev(- ex dz 


-f e 1 $ 1 есе) 
© J-a арт ау її 
@- х)? 
x exp(- p ) lz 
and a B E Ge a á 
da 000 ая ан)“ 


EL 


таи 
7а at 


(х, t) as required. 
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31.2.3 Summary of Section 31.2 


In this section we defined the terms 
repeated integral (page C30) 
convolution of two functions (раве C31) 
Theorem 
(Convolution Theorem, page C31) 
If 


Techniques 


1. Use the Fourier transform to solve partial differential equations in 
two variables, in particular the one-dimensional wave equation and 


heat conduction equations. 
2. Reverse repeated infinite integrals. 
Notation 


$. (page C25) 
y, (page C25) 
frg (page C30) 
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31.3 FOURIER TRANSFORMS AND EUCLIDEAN 
SPACES 


31.3.1 Fourier Transforms and Euclidean Spaces 


To round off this discussion we shall show how Fourier transforms are 
related to the Euclidean space theory we were studying when we first came 
across the idea of a Fourier coefficient (Unit 19, Least-squares Approxi- 
mation, pages K286—7). Since we have only considered real Euclidean 
spaces we shall confine the discussion mainly to the types of Fourier 
transforms which are guaranteed to be real, the sine and cosine transforms 
F, and F, defined in sub-section 31.1.1. To make the relation between a 
function and its transform symmetric we shall modify the factor 2/1 which 
appeared in the inversion formula for our previous definition. That is, we 
re-define the cosine transform of f to be 


5 5 fo 
Ao- E FQ) = J 2 | Го) cos (xy)dx (YeR) 


л х Јо 


so that the inversion formula is 


f(x)= К. F ЛО) соз (xy) dy (re [0, oo) 


т Jo 
A set of conditions under which this pair of formulas is valid is: (i) f has 
domain [0, оо) and codomain А; (ii) fis continuous and /' and /” are piece- 
а 
wise continuous; (її) f is absolutely integrable on [0, co), i.e. | ЛО dx 
о 


converges. The set of all functions with these properties forms a vector 
space which we shall call V[0, со). We can attempt to convert this into a 
Euclidean space in the same way as we did for the other vector spaces of 
functions we have dealt with, by defining an inner product on it. The 
obvious definition to use is 


© 
i-g= | p 
o 
with the associated Euclidean norm |, where 
E 
i-re | f*. 
o 
Before proceeding, we must ask whether this definition really does give 
us a Euclidean space. First of all, does f · g exist for every f, g in V[0, co)? 


ә 
The answer is, it does not: there are functions in V[0, оо) for which | f 
o 


does not exist. An example is given in Appendix 2. It follows that V[0, оо) 
with the above norm is not a Euclidean space. A similar situation arose in 
Unit 20, Euclidean Spaces II, where we wanted to define a Euclidean space 
structure on the vector space R” consisting of alhinfinite sequences of real 
numbers. We can deal with the problem in the same way as we did there: 
we define our Euclidean space to be the subspace of V[0, oo) for which 


re 


does converge. 


Let us call this space E[0, со). It can be verified that the axioms for a vector 
space are satisfied, so that E[0, со) is a vector subspace of V[0, со), and 
that f · g exists for all f, g in E[0, co). Thus E[0, œ) is a Euclidean space. 
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The question naturally arises now whether the Fourier cosine transform 
(i.e. the linear transformation of E[0, oo) to itself 

Si— fd 
is a Euclidean space morphism; that is, whether f: g = {с g, so that the 
diagram 


сабе 


NV 


fa 21с:9с 


commutes. Using the Fourier inversion formula we find, in fact, that 


tem [feet a 
" | Í * о) (fio cos (xy) dy} dx 
л Jo о 
2 | а 
= EI (f TOEO) cos (9) di] dy, 
n Jo о 


provided the conditions for reversing а repeated integral are satisfied. 
Subject to these conditions, we may therefore conclude, using the definition 
of f., that 


гв Í ЛОО) ау 


= К 
ѕо that 
7—7 
is a Euclidean space morphism. 


The formula we have just obtained may be written 
Frowa- |" zozoa 
or, in the special case where g =f 
| Yeo? dx= is LOM dy 0 


This last is called Parseval's formula, since it is an analogue of Parseval's 
equality given in Unit 20 (page K319) 


& 
lxi? = У (х e)?. 
06% 


We can derive Parseval's formula (1) from Parseval's equality for Fourier 
cosine series оп [0, р) as follows. From Corollary 9-1 on page K363, 
Parseval's equality for a Fourier cosine series takes the form 


Р 2 
f 0) ах =} (+a +ай+ ж) О] 
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where 
2 (» k 
«= | Г) cos == dx. 
р Јо р 


If we take the limit for large р of (2), an argument similar to that in sub- 
section 31.1.1 turns (2) into (1). 


Exercises 


l. Formulate the analogue for Fourier sine transforms of the Parseval 
formula for cosine transforms obtained above. 

2. Show (ie. with the same lack of rigour as in the above proof of 
Parseval's formula) that for real functions f, g with domain (— co, co) 


fg =f" JOO) dy 


= = 
where {+ g now means | S(x)g(x) dx, and g(y) means the complex 
E 


conjugate of g(y). 
3. We saw in the example of sub-section 31.1.3 that the Fourier trans- 
form of 
X — 25, etl jg y —_, E LT 
па +y? 
Use this fact, together with the version of Parseval's formula obtained 
in the preceding exercise, to evaluate 


= dy 
-o (a? + у*)ф* + у?) 
Solutions 


1 NIS deu 


where 


* ROO) dy, 


pom / 2 [лә sin (xy) dx 
Tt Jo 


can be obtained in the same way as the formula for cosine 
transforms. 

2. Since g is a real function, taking the complex conjugate of 
the result of applying the Fourier inversion formula 


ы” i; 
#0)е7 ° dy, 


1 
80) = —= 
Мт 


gives 


= 


= d = 
8(х) = g(x) = V | Bye" dy 


and so 


(assuming that the repeated integral may be reversed) 


f-g= Г fex [v dy) dx 


= ji JOZO) dy. 


35 


by the definition of f. (The last integral above is an example 
of a Hermitian form, as defined in Section IV-12 of N for 
complex vector spaces. Such forms give the generalization of 
Euclidean space theory to complex vector spaces, such as 
the one in which /'and g lie.) 

Writing f(x) for е7 1“! and g(x) for e^, we have 


» d 
-a (a? + у)? + y?) 
_ {7 EID FEN, 
i 520) ER) % 
ИЕ 
=й f TORO) dy 
af? 
-5f Sga) dx 


(by Parseval's formula, Exercise 2) 


xem. 
= zl есею dy 
2ab J-a 


а 
= — gets 


~ abla +b)" 


31.3.2 Summary of Section 31.3 


In this section we defined the terms 


Fourier cosine transform (re-scaled) (page C33) 
Parseval's formula (page C34) 


Technique 


Use Parseval's formula in the evaluation of certain integrals. 


Notation 
V[0, со) (page C33) 
Е(0, oo) (page C33) 


Main Result 


The mappings f — J, and f———— j, with domain £[0,co) are 


Euclidean space morphisms.” 
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31.4 SUMMARY OF THE UNIT 


This unit introduced: another technique for solving partial differential 
equations: the Fourier transform method. There are two approaches to 
the definition of a Fourier transform. One is through investigating the 
Laplace transform of complex functions and the other is by investigating 
what happens to the convergence properties of Fourier series when the 
size of the interval on which they are defined is extended. It was this latter 
approach which we used in the unit, although we did subsequently inves- 
tigate the connection between Fourier and Laplace transforms. In this 
way we defined the Fourier sine and cosine transforms, and these led to 
the definition of the Fourier transform f, where 


fu) = X Í Ёз Sle" ах. 


The integrand in this definition is a function from the reals to the complex 
numbers. We therefore investigated various techniques of integration for 
such functions. 


In the second section we illustrated the properties of the Fourier transform 
by using it to solve the general one-dimensional wave equation and the 
heat conduction equation. Since the solution to such equations involves 
integrals of the form 


EI 
| ЈО) x Ben" dy 
Е 
where f and g are known, we were led to define the convolution of fand g 
denoted by f » g. 
1 * 1 = 
feste T [7 ree nm |? лада 
Jin AES 2n J-« 
The convolution theorem states: 
if 
h-fxg 
then 
h=frg. 


In the final section we set up the vector space V[0, со) of all functions 
which 


(i) have domain [0, со) and codomain R, 
(i) аге continuous and have piecewise continuous first and second 
derived functions, 
(iii) are absolutely integrable on [0, со). 


In И[0, со) 

Јо) = E [ (х) cos (ху) х — (ye А) 
апа ‚ 

Јо) = E [ ЛО) соз (ху) dy (хе [0, о) 
The subspace £[0, co) of V[0, oo) for which 


l^ 
o 


converges is a Euclidean space with inner product 


га | 
b 
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The linear transformation 
f——f. (fe El, 0) 

is a Euclidean space morphism, i.e. 
Pesk è 

Definitions 
finitary function (page C6) 
Fourier cosine transform (pages C7 and C33) 
inversion formula (page C8) 
Fourier sine transform (page C8) 
integral transform (page C10) 
Fourier transform (page C16) * 
absolutely integrable function (page C16) 
repeated integral (page C30) 
convolution of two functions (раве C31) 
Parseval's formula (page C34) 


Theorem 


(Convolution Theorem; page C31) 


If 

h-fxg 
then 

h=feg 
where 

fas |7 Лднк-э& 

1 fe 
= Tau [Г = z)g(z) dz. 

Techniques 


l. Integrate functions from A to C by applying the rule 


b 'b b b 
[4] iit = [ veil k 


where g = h + ik. 
2. Apply the following techniques of integration to functions from R to C: 
(i) integration by substitution, 
(ii) integration by parts, 
(iii) differentiation under the integral sign. 

3. Use the Fourier transform to solve partial differential equations in 
two variables, in particular the one-dimensional wave equation and 
heat conduction equations. 

4. Reverse repeated infinite integrals. 

5. Use Parseval's formula in the evaluation of certain integrals. 


Notation 
F, (page C7) 
Е (page C8) 
7 (page C16) 
Ф, (раве С25) 
V, (page C25) 
f*g (page C30) 
У[0, оо) (раве C33) 
E[0,co) (раве C33) 
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31.5 SELF-ASSESSMENT 


Self-assessment Test 


This Self-assessment Test is designed to help you test your understanding 
of the unit. It can also be used, together with the summary of the unit, 
for revision. The answers to these questions will be found on the next non- 


facing page. We Suggest that you complete the whole test before looking 
at the answers. 


І. Find the Fourier transform of f where 
-fb if |x| <a 
fe = (6 if [|х|>а] CER 
2. (i) Assuming that the Fourier inversion formula applies to the 
function in Question 1 despite the fact that f is discontinuous 
at +a, evaluate 
* si y 
| in (ау) cos (ух) dy. 
-=o y 
for |x| # a. 
(1) What value would you expect for the integral in the case x = а? 
* si 2 
3. Evaluate Í m Кл dy. 
e y 
(Hint: Use Question 1 and Parseval’s formula for Fourier integrals: 


[vera [* var 


4. Evaluate: 


н a. 
() dy f sin (xy) dy, 


b 
(ii) Í уе? dy. 


5. Explain the rôle that convolution could play in the evaluation of an 
integral in the form 


19)= | SO) x g0)e™™ dy 


6. The differential equation that models transverse vibrations of an 


elastic rod is 
д%и 1 ĝu 
D 


where и: R? ——+ R. 


Gi) Without justifying the steps in your working rigorously, use the 
method of Fourier transforms to obtain a general solution of 
this differential equation (analogous to the solution of the wave 
equation in Equation (7) on page C27) as the sum of two integrals 
each involving one arbitrary function with domain R and co- 
domain C (the field of all complex numbers). 

Gi) In the case where П 


u(x, 0) = A(x) 
ди (хє А) 
a 9-0 


with ^ a known real function, deduce that a solution to the 
problem is 


u(x, D- Z [| F [ |. Axel? as] cos (ay^tje^* dy. 
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Solutions to Self-assessment Test 
1. The Fourier transform is f where 


= = | Slade”? dx 
2л 


belt” I 
= | -= (y #0) 
рст -a 
b fee 
( iy 
с. 2b e ay! 
Е y 


) (у # 0). 


1 fe 1 fe 
= dx = bdx= =, 
700) Fal fe A]. dx Ta 


2. (i) The Fourier inversion formula is 
1 (= zi 
ЈО) = xl Fe" ay. 
Ул -o 
From Question 1 


1 = 2_ (2) а= 1 |x| «a 
Mn] -= J2n V y о |х| >а 


тее бею), 
т Ј-ә y 


zi sin (ay) 809 ay - | sin (ay) sin (xy) dj. 


RJ-o P, aj- » 

The integrand in the second integral is odd and so the integral 
is zero. 

Hence 


© sinaycosyx | [n |x| «a 
Cee eG Bis 
(ii) At the point of discontinuity x = a, we expect, by analogy with 
Theorem 9-1 on page K340, the Fourier inversion formula to 
converge to the mean of f(a*) and f(a~), which is jb. 
Putting b = 1 we therefore expect 
l[* sin ау соѕау 
Hi ——— —— dy 5d. 
jas | У 
Hence ` 
|с sin ay cos ay 
=o » 
3. From Question 1 with 5 = 1, 


jo- 2 e 2) 


NIN 


Jan\ > 
when fo) = (о oe 


Using Parseval's formula this gives 


а -2: fsinZiay 
[rove [IE 


LM 31.5 


But 


[ (1)? dx = 2а 


b 
= | x cos (xy) dy 


= [sin (xy) 
= sin (xb) — sin (xa). 


M b d [е 
ii 7 IXY d) - =j 
Gi) [= у [>> (©) y 
pel] b pixy 
= p ] -Í dy 
IX Ја а IX 
(by integration by parts) 


=) eit, 


If f and g can be recognized as the Fourier transforms of two known 
functions, л and К say, then 


LE к, f 
Ux) = Í Щу) x k(y)e7^? dy. 
-® 
We can then apply the convolution theorem to obtain 


нй йз Г h(z)k(x — 2) dz 
Уа) 


zc p щх — 2)k(2) dz. 


Эл 


(i) The transformed equation is 
1 


Jon 


o gtu | 
| эя ® Deh” dy 


-% 
3 zh 2 Е (х, Ne" dy. 
If we introduce functions y;, (one for each y) and ф, (one for 
each г) where 

$, =x ——9 u(x, t) (хє А) 
and 


dy:y—— G0) WER) 
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the transformed equation becomes (using integration by parts four 
times on the left-hand side and differentiation under the integral 
on the right-hand side) 


(ayy = = 3 9500 


Wy = = (ay? PY). | 
This ordinary differential equation has general solution (one for 
each y) 


ry) = Ge"! + Рет", 

We generalize the arbitrary constants С, and F, by'writing 
Vt) = ен" + /o)e7 

and then use 


VO = à) 


and 


uie —— р e^ "G (y) dy 


Ул 


to obtain 


ue, =e [f етее у) dy 
J2n Ј-ә 


+ E g tom) «| 


T es à MEC 
(ii) Substituting for u(x, 0) and > (x,0) (using integration under the 


integral sign in the second case) we find 


1 е Ld È 
Meh = = [ | Eae [^ Foye" dy] 


1 B è 
0=—— [| #0)іау?ет 9 dy 
V2n Ш-ә 
EI 
+ | лао] 
-9 
If f and g are the inverse Fourier transforms of Гапа g, then the first 
formula gives 
A) =g) + О) 
and the second can be satisfied by making 


ЈО) = #0) 


Јо) = (ә). 
Hence 
f-godh 


Substituting this in the solution obtained in (i) and using Euler's 
formula 


(е0 + е7) = cos 0 


with 0 = ay*t, we obtain the given formula. 
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316 APPENDICES (Optional) 
Appendix 1: proof of formula (7) in sub-section 31.1.1. 


(We expect this to be of interest primarily to those students who enjoy 
analysis.) 


We wish to prove that, for any x>0, 


T g= lita h(g(a) + g(a + h) + g(a + 2) +} a) 


where 


80) = cos (xy) |. 2) cos (zy) dz o 


We have assumed that f is continuous, piecewise twice continuously 
differentiable, and finitary. Hence in [0, со) /” is continuous except at a 
finite set of points, which we call, in increasing order, 


Xp Xas eses Xag 
where 

хо = 0, x, = Po 
with [— po, po] the interval outside which f(x) = 0. 
We shall prove three lemmas 


b 
(i) g is continuous (this ensures шш [ g exists). 
a 


(ii) There is a number A such that 
А 
0) < = 
Е] m 


b 
for all y (this ensures that lim | g exists). 
boo Ја 


(їй) The “tail” of the series, i.c. 
T, = h(g(a + nh) + g(a + (n + 1)h) + ses. 


is small for large n (this ensures convergence of the series on the 
right of (1)). 


Proof of (i) 

is the product of the continuous function 
y (———35 cos (xy) 

and the function 


F.:ig——— ju J (x) cos (xy) dx, 
° 


so it is sufficient to prove that F, is continuous. This is done as follows: 
for any є we have 


FO +0) = FQ) = [ros x(y + 0) — cos (xy)] dx 
° 


2, -2 [ro sin x(y + 3c) sin (4xe) dx 
ó 


since 
2 sin г sin 5 = cos (r — 5) — cos (r +5) 


for all real r, s. 
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Непсе 
IF. )- 501 2 firo [sin хоу + 4e) [sin xo] d 
o 
n 
<2 f |JGO|ixedx since |sin 0] <19] 
Јо 
ro 
< poe [st] s 
" 
and hence 


lim Ку + €) = Fo) 
e~o 


for all y, which is the definition of continuity. 


Proof of (ii) 
Integration by parts gives 


к) = fo cos (ху) dx 
e» | f) cos (ху) dx 
i sos 


ЕС ү sin (xy) „n Sin (xy) 
= [Œ y - fre y a] 


mex Í = SW sin (xy) dx 
Y= Jains 


En 


since sin 0 = 0 and f(x,) = 0 and the other terms not involving integrals 
cancel in pairs because of the continuity of f. 


A second integration by parts gives 


1 


n * 
ЕДу) = – 7 à [7 '(X) eos (xy) + I “(x) cos (xy) а] 
i= En 
This time there is no cancellation since f* need not be continuous at 
Хү, Xz) +++) Ж-ү, but the magnitude of the ith square bracket is bounded 
above by 


xi 
Ire cerea [^ urere 
xii 
and therefore the whole sum is bounded in magnitude by the sum of 
these expressions, which we may call 4. We conclude, using (2), that 
A 
ID <| FO sur 


Proof of (iii) 
17, «Аба + п) + |ба + (n + 0) +) 
А А а 

“com ерту) "o 

< їй h A dy (see diagram) 
Ja+(n—1) Y 

a A 

с а+(л—1)/ 


where we have assumed that n is large enough to make a+ (n — 1)һ> 0. 


The area of the shaded blocks 
represents the sum, which is 
less than the area under the 
Curve, representing the integral. 


а+{п-1)һ anh а+(п+1)Һ y 


Using these three lemmas we can now prove (1) itself. We shall show that, 
for large b and n 


[р gx [ emm nnm entm 1) 


ж h(g(a) + gla +) +) [9] 


where = (b — a)/n and thesymbol ~ means that the expressions it connects 
can be made as close together as we please by choosing b and n large 
enough. 


From Lemma (ii) we have 


а b E © СУ А 
Fee [fe ma faje [Saad 
a a b b b b 


Also, by lemma (iii) we know that T,, the difference between the finite 
and the infinite series in (3), satisfies 


A A. Hk 
ёс зы e M ud 
тета вт ti 


Thus, by choosing b large enough, we can make both of these differences 
as small as we please provided that / is less than 1. The remaining difference 
to consider is 


b 
Í g — h(g(a) +++: + g(a + (n — 1)h)) 


and, once b has been chosen, we can make this difference as small as 
b 
we please by making n large enough. (This is just the definition af ГА 


given іп Unit M100 9, Integration I.) Thus by suitably choosing first b and 
then п we can make all the members of (3) as close together as we please. 
Since the first and last members do not depend on 5 огл, this is only 
possible if these two members are equal. 
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Appendix 2 


A function that is absolutely integrable, but not square-integrable, on 
[0, оо) is 


where 


В 1 
inet 
) 

\ 


(x e [0, оо)) 
1-4 ба - n) ifn ёки; 
т л 55646 25 
1 fix) 
т > 
1 2 3x 


The graph of f consists of a succession of parabolic arches, the nth of 
F А KRR s д " 
which has height n and width 5, so that its area is proportional to 
п 


1 2 
arm actually <3 s 


Thus we have 
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The graph of f?, on the other hand, consists of a succession of arches of 
x es É " 
height л? and width 73; the integral 
n 


MEE ag 
i". 


is therefore proportional to 
x 1 A 


© | а ә 
and since У, ~ diverges the integral | J? diverges, and hence | ў 
п + 


п=1 o 


diverges too. 
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